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1 Introduction 



The final state of a heavy ion colhsion has been recently studied in gauge/gravity duality in a 
number of papers P [3 (3 ll El [71 El [9l [TOl [III [IS [II El [l5] . The early stages of the collision 
process are more difficult to model and have so far received less attention [16]. In this paper, 
we consider both the final state and the onset of the collision in a simplified model in which the 
boundary theory is (1+1) dimensional and the bulk is locally AdSs. We shall give explicitly 
the general bulk solution that is dual to two colliding clusters of conformal boundary matter. 
That such a general solution can be found is, of course, due to the simplicity of gravitational 
dynamics in three dimensions. The solution involves two arbitrary functions, each depending 
on only one of the light cone coordinates. With a suitable choice for these two functions, 
the boundary process can be described in the language of heavy ion collisions as having two 
incident extended systems, an interaction diamond and central and fragmentation regions in 
the final state. There is one important difference between our boundary matter and that 
of real (3+l)-dimensional heavy ion collisions though: our matter is extremely opaque and 
bounces off backwards in the collision. 

As a special case, our bulk solution also provides the metric of two colliding shock waves 
in three dimensions. This metric describes a spacetime in which two massless point particles 
create a three-dimensional black hole \n\ [T8] . In four (or more) spacetime dimensions the 
metric of a single shock wave is known both with a vanishing [191 [20] and nonvanishing 
[2 H 1221 [23] cosmological constant, and an extension to colliding waves would have important 
applications in the study of black hole production at colliders [23] , but finding such a solution 
has proved a formidable task j25[ [26] . 

We use the three-dimensional coordinates {x^,x~,z), {T,r],z) and {t,x,z), referred to 
respectively as the light cone, Milne and Minkowski coordinates. The conformal boundary is 
at 2 = 0, the coordinates are at z = related to each other by 

x^ = — ^ = —;=e^'^' t = T coshri, X = T sinhn, (1) 

where r > 0, and the boundary Minkowski metric takes the form 

ds^ = -2dx+dx' = -dr^ + r^dry^ = -dt^ + dx^ . (2) 

2 Colliding conformal boundary matter 

In three dimensions, the metric of a plane wave moving in the direction can be written 
in the light cone coordinates as 



9mn — ^2 



/ -1 

-1 zV(x-) 1, (3) 
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where f{x~) is an arbitrary function of x~ . This metric satisfies the AdSs Einstein equations 

Rmn — \ RdMN — -^gMN = 0, (4) 
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following from the action 



1 r 3 . 2\ 



d^xV=^LR+— , (5) 



9mn — ^2 



and must therefore be locally isometric to AdSs. In the special case in which f{x~) is propor- 
tional to 5{x~)^ we obtain the three-dimensional negative cosmological constant generalisation 
of the Aichelburg-Sexl (AS) shock wave |19j . 

By symmetry, the metric of a plane wave moving in the x~ direction can be written as 

/ z^g{x+) -1 \ 

-10 0, (6) 

V 1/ 

where g{x^) is an arbitrary function. Note that in each case the metric components are 
independent of the light cone coordinate in whose direction the wave is travelling. 

While extensions of the above plane wave solutions to four and more dimensions are known 
|2 H I22 t [25], an extension that would describe two colliding waves in four or more dimensions 
has proved elusive. In three dimensions, however, we can find such a solution by making in 
the light cone coordinates the ansatz 



9MN = -T r 1 h (7) 




where the two-dimensional metric components g^j^y may a priori depend on any of the coordi- 
nates. If g^y is assumed to reduce at z — > to the (l+l)-dimensional Minkowski metric ([2]), 
the general solution to Einstein's equations @ may be shown to be 

2 / z^9{x+) -1 - 

g{x+)f{x~) z'f{x~) I , (8) 



9MN = 
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where g and / are arbitrary functions, each of dimension 1/ (distance)^. This solution can also 
be obtained by specialising that given in [27j to a flat boundary metric. It may be regarded 
as a superposition of the plane wave ([3]) moving in the x~^ direction and the plane wave Q 
moving in the x~ direction, the only nonlinearity in the superposition being a correction in 

the component g^ The solution can thus be characterised as a collision of two plane waves. 

If / and g are each proportional to the delta-function, the solution describes the collision of 
two shock waves. 

To analyze the bulk solution ([8]) from the boundary theory viewpoint, we compute the 
energy-momentum tensor on the conformal boundary at z = [28]. The small z expansion 
of g^j^u reads 



g^y{x^,z)=g^l^} +g^^}z'^ + ... 




\9{x+)f{x , , . 

\g{x+)f{x~) 



(9) 
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from which we obtain 



T^jy is clearly conserved. A conserved total energy E can hence be defined as the integral of 
T** over a surface of constant t. From (jlOp we find that E decomposes as 



oo 



E = E+ + E-, 
E-- = —pT^ / dx+g{x+), (11) 



where and E~ are respectively the energies of the components moving in the direction of 

x'^ and x^ . 

When 5 = but / 7^ 0, the energy-momentum tensor can be written in terms of the null 
vector = (1,0) as 

which shows that the boundary matter is null dust moving in the direction. Similarly, 
when / = but g ^ 0, the boundary matter is null dust moving in the x~ direction. When 
both / and g are nonvanishing, the energy-momentum tensor is the sum of the left-moving and 
right-moving parts, and as / and g are independent, it is possible to interpret the boundary 
matter as two null dust clouds that just pass through each other without interacting. In 
higher dimensions interactions are however present, and the boundary energy-momentum 
tensor after the collision is then usually interpreted as a single-component fluid with internal 
dynamics and thermodynamics [29]. As we are viewing the (1 -|- l)-dimensional boundary as 
a simplified model for the (3 -|- l)-dimensional boundary, we shall seek a single-component 
interpretation for the boundary matter. 

Consider in particular the situation in which / and g are both non-zero and have the 
same sign. We can then write the energy-momentum tensor in the perfect fluid form T^,y = 
(e + p)u^Uy + pgliu , where the normalised future-pointing timelike velocity vector reads 



g{x+) f f{x ' ^ "^^^^ 



\4f{x-) J ' \4g{x^ 
and the energy density e and the pressure p are given by 



(13) 



P 



'8ttG? 



^g{x+)f{x-), (14) 



with the sign in (jl4p being that of / and g. Note that the dimensions in these formulas are 
correct since C/G3 is dimensionless and / and g each have dimension l/(distance)^. To find 
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the trajectories of the comoving fluid elements, we have to integrate the equation dx^ /dr = 
u^{x^,x^). Using the normalisation condition 1 = —u^ = 2u^u^ = 2{dx^ /dT)'^{dx^ /dx^), 
we find = (^^dx~^ /2dx~ , ^/dx^~/2dx^^ , and comparison with (jl3p then shows that the 
trajectories are obtained by integrating the separable differential equation 



dx _ / ff(x+) 

d^-^Ji^)- ^^^^ 

3 Special cases 

While the general bulk solution ^ involves the two arbitrary functions / and g, Einstein's 
equations ^ imply that this solution must be locally isometric to AdSa. The solution 
can therefore be brought at least locally to the standard Bahados-Teitelboim-Zanelli (BTZ) 
form [301 EU, 

where 



ds^ = -Fdr + — + r\dip + N'^dt), (16) 
F 



F= ^-'-^y-\ «' = -^. (17) 

and the the parameters r± determine the mass parameter M and the angular momentum 
parameter J by 

We shall now write out this transformation for certain choices for / and g. These choices 
will emerge in section H] as regions of special interest in a solution that describes a heavy ion 
collision process on the boundary. 

3.1 Bjorken similarity flow 

Consider the region in which > 0, and let 

g{x) = fix) = , (19) 

where M is a constant. For M 7^ 1, we have from (jl3p 




= \-\li::^^-\li:::r]=^i-e-^,-e% = ^ = i=(e^e-''), (20) 



and (IHI) gives 

_ C M-1 _ C M -1 

In the coordinates (t, x), (j20p reads = {I/t^x/t). The boundary flow is thus Bjorken 
similarity flow [32J. 
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In Milne coordinates {T,rj,z), the bulk metric reads 

22 



ds = — 



£2 

y2 



{M-l)z^\\ 2 , A , (M-l)z2\2 



This bulk metric was analysed from the gauge/gravity duality viewpoint in [2]. The metric 
can be transformed to the BTZ form (jl6p with J = 0, 

^-Mjdt2 + -3^^^+rW. (23) 

For positive M, the Killing vector 9f has a Killing horizon at r = L\fM^ and the temperature 
associated with this horizon is T = \fM j i^ixC) and the Bekenstein-Hawking entropy A^btz 
for an interval At/ is given by A5BTz/(£^f?) = \/M/(4G3). A time-dependent scaling 
argument from the boundary of (j23p to the boundary of (j22p then gives the flow (j20p - (j2ip 
the time-dependent temperature and entropy density 

= — , = • (24) 

27rT 4G3 r 

The term —1 in the factor M — 1 in (|2ip can be interpreted as a vacuum energy contribution, 
and once this contribution is subtracted, the fluid is a perfect fluid in adiabatic expansion for 
any positive value of M. 

3.2 Angular momentum 

Let /(x) and g(x) be again proportional to 1/^2, but with coefficients that are not necessarily 
equal. We write the coefficients as 

4(x+) 4(x-) 

where M and J are constants. In the domain > 0, going to Milne coordinates (r, rj, z) then 
puts the metric in the form investigated in section 5 of [9j, and it follows from the analysis 
therein that the metric can be brought to the spinning BTZ form (|16p . If J/ C = M — 1, so 
that g vanishes and we have on the boundary a null dust p2p moving in the direction, 
this conclusion can be extended from the domain > to the half-space x~ > 0. Similarly, 
if —J/C = M — 1, so that / vanishes and we have on the boundary a null dust p2p moving 
in the x~ direction, this conclusion holds in the half-space > 0. 

3.3 Plateau: constant / and g 

Let now both / and g be constants. 

Suppose first that / and g either have the same sign or are both zero. By a boost in 
(x+, x~) it is then always possible to make / and g equal, so it suffices to consider this case. 
We set f = g = fi/{2C'^), where is a dimensionless constant. The bulk metric reads 



ds^ 




{dx+y + (dx-) +dz'} , (26) 
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and transforming to the coordinates {t,x,z), brings this metric to the form 



AC 



dr + C^^+[- + ^] dx\ 



C ^z\' 



AC 



(27) 



Writing finally 



C 



C fiz 
l^AC 



the metric is transformed to the spinless BTZ form (|23p with dr] 
/i > 0, the temperature associated with the Killing vector dt is 

27r£' 

and the entropy ASbtz for an interval Ax is 

AS'btz _ ^/M 



Ax 



AGs 



(28) 

dx/C and M = fi. If 
(29) 

(30) 



If / and g have opposite sign, a similar analysis can be given, but this case will not occur 
in section H] below. 

Suppose then that g = but / 7^ 0, so that on the boundary we have the null dust (|12p 
moving in the x"*" direction. We now write / = [ij C? , where /Li is a nonvanishing constant. 
The bulk metric is 



£2 

ds^ = 

Z'^ 



-2dx^ dx + ^-rrridx f' + dz^ 



£1 ^ 

z2 + 2 



dx 



cyz^ + fi/2j 



(31) 



and can, by introducing the coordinates {t, r, x) with 



be written in the form 



ds^ 



(r2-^£2/2)' 



dx {Cfi/2)dt\' 

z ^ ) 



(32) 



(33) 



This metric is of the BTZ form (fTHl) with dip dx/C, M = fi and J = For /i > the 
metric is hence an extremal rotating BTZ hole, with the angular dimension unwrapped. 

A similar analysis applies to the case in which / = but g ^ 0, with M and J now having 
opposite signs. 
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3.4 Interpolation 



Finally, we set 

fix-) = = ^7^, (34) 

where a is a positive constant of dimension length, t> is a positive dimensionless constant 
and we assume x"*" > 0. This case can be regarded as being half-way between the Bjorken 
similarity flow case (I19p and the plateau case of subsection 13.31 Note that d^- is a spacelike 
Killing vector in the bulk, and the boundary flow (jl3p is invariant under the corresponding 
null Killing vector d^- on the boundary. 

We first transform from {x~^,x~,z) to {w~^,w'',a) by 



X 



vaw 



log 



log ( J 



1 + (vl2) 



: W 



-a + 



{vl2) 



1 + (yl2) 



: W 



(35) 



The conformal boundary is then at a ^ co with fixed w'^ . These coordinates are adapted to 
the commuting bulk Killing vectors = av d^- and 9^+ . To eliminate from the metric a 
cross term proportional to dadw~^ , we transform to {y^,y~,p), where 



cosh 2a = pJl + {v/2), 



dw 



dw 



dy' 
dy" 



{v/^)dp 



(p2_l)^/^l + (^/2)2 

(i>/4) pdp 



(p2-l) / l + (^/2)^ p^-l 



(36) 



and p > 1. The metric becomes 
ds^ = 



{dy+f + {dy-f-2pdy+dy-+ '^^ 



A{p^-l) 



(37) 



Note that both a and v have now disappeared. Transforming further to {9^ , 0^, x) by 



y 



cosh 2x, 



where x > 0, we have 

ds^ = C? 



- (sinh xf{dB^f^ (cosh x) ^ (d^^ )% d^ 



(38) 



(39) 
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The metric (j39p is recognised as part of the AdSs hyperboloid in coordinates in which the 
Kilhng vectors dgi and 5^2 generate two commuting boosts [30l[3l]: writing the hyperboloid 
in embedding coordinates (T^, T^, X^, X^) as 



,2 / „\ 2 / , \ 2 



ds^ = -(dT^) - (dT^) +[dX^^ + (dX^) , (40) 

we can choose the transformation to (f39l) so that dgi = T^dx^ + X^drpi and 9512 = T'^dx^ + 
X'^dx^- As d^- = {va)~^dy~ = {va)~^ {dgi — dg2), it finally follows that 

d^- = — (r^dxi + X^dri - T^dx2 - X^dT2) . (41) 



va 

To summarise, the invariant geometric characterisation of the bulk Killing vector d^- can 
be read off from (j41|) . As d^- is the bulk continuation of the boundary Killing vector under 
which the fluid flow is invariant, one is tempted to seek thermodynamics for the boundary 
flow from thermodynamics of bulk phenomena that are invariant under d^- ■ This strategy 
works for Bjorken similarity flow, as outlined in subsection 13.11 [9], and for the plateau flow, 
as outlined in subsection 13.31 since in both cases the Killing vector in question is a spacelike 
boost and the Killing horizon of the commuting orthogonal boost has a Hawking temperature. 
With the Killing vector (14ip . however, a thermodynamical bulk interpretation appears not 
to be known, and there is evidence that one may not exist |33| . It is certainly possible to 
transform the bulk metric to the BTZ form: given any M and J satisfying M > \ J/C\, the 
transformation 

1 



2/± = -^M^{J/C)^-±ip 
2{r/Cf - M 



'M^ - {j/cy 

takes the metric (1371) to the BTZ metric (1161). However, we then have 



(42) 



Cdt±d^ 



and consideration of the Killing vector d^- = {va)~^dy- does not appear to single out specific 
values of M and J for which the Killing horizon of dt could be argued to provide a temperature 
for the boundary flow. 



4 Collision 

We now turn to situations in which a collision of two null dust flows evolves into Bjorken 
similarity flow at late times. We hence assume both / and g to vanish at negative argument, 
be nonzero and have the same sign for positive argument and approach Bjorken similarity 
flow form ()19p at large positive argument. 
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Figure 1: Structure of the conformally invariant boundary flow following from the solution 
assuming that / and g are vanishing at negative argument and are nonvanishing and have the same 
sign for positive argument. Within the past light cone of the origin the energy- momentum tensor 
vanishes. The left wedge < < x~ contains null dust moving to the right, and the right wedge 
< < contains null dust moving to the left. The future light cone of the origin contains 
expanding matter with the flow given by (jl3p and energy density given by (jl4[) . The special solution 
(|44ll that gives rise to unregulated Bjorken similarity flow ([20| - ((2T|) within the future light cone of 
the origin is also indicated; for this solution the energy density ([2T|) diverges as r — > 0. 



The general structure of the boundary flow is shown in Figure [TJ Within the past light 
cone of the origin we have vacuum, the left wedge < < x~ contains null dust moving 
to the right, and the right wedge x~ < < x^ contains null dust moving to the left. Within 
the future light cone of the origin the null dusts have collided and formed the flow given by 
(jl3[) and (I14p . which at late times tends to Bjorken similarity flow, ()20p and (I2ip . 

The key point is now that a suitable choice of the null dust profiles resolves the initial 
singularity of Bjorken similarity flow. We look at unregulated Bjorken similarity flow in 
subsection 14.11 analyse two regulated versions in subsections 14.21 and 14.31 



4.1 Unregulated Bjorken similarity flow 

Consider first the choice 

9{x) = f{x) = ^^e{x), (44) 

where M is a constant. Within the future light cone of the origin we then have Bjorken 
similarity flow ()20p - (j2ip . As reviewed in subsection 13. 1^ for positive M this flow can be 
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understood to describe a perfect fluid in adiabatic expansion with e(r) ~ and T{t) ~ 
1/r. 

The problem with this scenario is however that the energy density diverges at r — > 0, and 
the energies (llip of the initial null dust components also diverge. While for a (3+1)- 
dimensional boundary the bulk dynamics may provide criteria for regulating the initial sin- 
gularity [3 [13], we are not aware of similar criteria for (l+l)-dimensional boundary. We shall 
therefore regulate the singularity by a bottom-up phenomenological approach, introducing in 
the functions / and g and new scale that corresponds to the initial conditions. 



4.2 Smooth profile 

As a first alternative, we modify Bjorken similarity flow functions (I44p to 

/(x) = g(x)= 0(x), (45) 

where a is a positive constant of dimension length. As the only other length scale in the 
problem is £, one may expect on phenomenological grounds a to be proportional to >C by a 
dimensionless constant of order unity. For > 0, we then obtain from ^ 



e = p= ^ , = , ' (46) 

^ 8ttG3 v/(x+^ + a^)(x-^ + a^) lOvrGs ^t^ + 2a^ cosh{2r]) + 4a^ ^ ' 

The energy density thus remains finite at r — > 0, has some r] dependence at intermediate r 
and approaches the similarity solution (j21l) as r — > oo at fixed rj. 

To see how the singularity at r ^ is opened up, we may compute the flow lines. Now all 
the particles do not emanate from the origin x^ = x~ = 0. When the initial point is on the 
right-hand-side branch of the light cone, at = and x'^ = Xi > 0, equation (jlSp implies 



dx dx 



^ + a2 \ (x+f + a2 



(47) 



which can be integrated to 



x-^ + \l [x+f + a? ] - {xi + x xl: + a? 



x-{x+) = ^ ^ . > { • (48) 

2 {x+ + ^{x+f + a2 j fx, + y^xf + a2 j 

Note that these paths stay in the region x^ > x~ . When the initial point is on the left-hand- 
side branch of the light cone, the situation is similar with x"*" and x^ interchanged. The only 
path that starts from the origin is x"*" = x~, obtained from (j48p with Xj = 0. A selection of 
flow lines is shown in Figure [2j In the limit a — > 0, we obtain Bjorken similarity flow (j20p . 

What this regularisation seems to lack, however, is a way to associate a temperature to 
the early stages of the flow. With Bjorken similarity flow, the temperature arose by noting 
that the flow is invariant under longitudinal boosts generated by the Killing vector dr^, the 



10 





Figure 2: Flow lines ([48]) for a = 1 and Xi ~ 1/4, 1/3, 1/2, 1, 2, 3, 4 and 5, and the corresponding 
flow lines with x"*" and x" interchanged. 



bulk solution (j22p extends this Killing vector from the boundary to the bulk, and in the bulk 
metric (I23p there is a commuting timelike Killing vector dt whose horizon has a Hawking 
temperature [9j. We have not found a similar isometry argument that would apply to the 
flow with the particle paths ()48p . 



4.3 Hard core profile 

As a second alternative, we modify Bjorken similarity flow functions (j44p to 



/(x) = g{x) 



M - 1 



-^e(x - a) + ^0(x)6(a - x) 



(49) 



where a is again a positive constant of dimension length. Phenomenologically, we might call 
the second term in (09]) a hard core "nucleus", of characteristic size a, and the first term its 
longitudinal "tail". We shall assume M > 1. 
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The different regions on the boundary are shown in Figure [3l From (llip we find that the 
"incident beam energies" in the left and right wedges are given by 

and for each null dust component this energy is distributed equally between the nucleus and 
its tail. As shown in subsection 13.31 the bulk metric in each null dust nucleus region can 
be thought of an extremal BTZ black hole, with the angular dimension unwrapped, and has 
therefore vanishing temperature. We thus regard the incident null dust beams as completely 
unthermalised, both in the nucleus and in the tail. 

The region > 0, in which the beams have collided, is divided into the central region in 
which the null dust tails overlap, two fragmentation regions in which one nucleus overlaps 
the other's tail, and the interaction diamond in which the two nuclei overlap. In the light 
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cone coordinates, the velocity field is given by 



2x- 



2x+ 




2a ' V 2x+ 



2x- 




1 1 



— , central; 
r 



fragmentation at x < a; 



fragmentation at x"*" < a; 



diamond. 



(51) 



Note that in the interaction diamond the fluid is thus at rest in the Minkowski coordinates 
{t,x). For the energy density, (fill) and ([49l) give 



e{x~^,x ) 



£(M - 1) 
327rG3 



1 



-e(x+ -a)e(x- - a) 



H -G(x+ - a)e(x")e(a - x") + (x+ ^ x" 

ax+ 

+4T0(a - x+)e(x+)e(x~)e(a - X- 



(52) 



or, in Milne coordinates, 
£(M- 1) f 



??) 



327rG3 



e(T - aV2 



-^efiog ^ 



^/2 



H e \r]\ - log 



+e(a\/2 - 



1^/, aV2 
log 



+ Xf^e |r?|-log 



a^/2 
aV2^ 



ar 



T 



(53) 



In Eq. (j52p the first line describes the central region, the second line the two fragmentation 
regions and the last line the interaction diamond. In Eq. (j53p . plotted in Figure HI the first 
line gives the energy density for r > a\/2, where the interaction diamond does not contribute, 
and the second line gives the energy density for r < a\/2, where the central region does not 
contribute. At = 0, (j53]) reduces to 



e(T,r/ = 0) 



£(M- 1) 
327rG3 



2e(r-ay2) e(a^/2-T 



(54) 



We may again obtain a concrete view of how the singularity is opened up by computing 
the flow lines. For concreteness, consider the flow lines approaching from the right. In the 
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Figure 4: The energy density ([53]) divided by C{M — l)/(327rG3) as a function of t] at constant r in 
the regularisation (|49)) . The left panel shows the case r > a\/2 and the right panel the case r < a\f2. 
The parameter h equals a\/2. 

null dust region these lines are the null lines = xf, where xf is a positive constant. The 
lines continue into the other regions for j;^ < o as 



and for xf > a as 



X — , or X = 2;^/v2, diamond; 
x+ 

a — xf + a log — , fragmentation region; (55) 

x+e^i 1°- ^ central region, 



X 

alog— r, fragmentation region; 



X 



X- 

a 



(56) 

central region. 

ex- 

The lines are plotted in Figure [5l What is most striking here is that the fluid is brought 
momentarily to rest on the light cone of the origin, and in the fragmentation region it then 
accelerates in the direction opposite to where it came from. The conformal matter is thus 
extremely opaque. This seems superficially to be in complete disagreement with what one 
knows about heavy ion collisions, where nuclei are extremely transparent. However, this 
transparency is due to the valence quarks which represent non-conformal degrees of freedom 
in nuclei. By contrast, the conformal degrees of freedom are those built in the classical gluon 
fields (see Section [5]) computed from the valence currents as sources [34] . 

For further insight into the dynamics, we may consider the conserved energy (jlip on a 
surface of constant t. For t < 0, all of the energy is in the incident null dust flows, while at 
t = o/\/2 exactly half of the energy is in the interaction diamond and the other half is still 
in the null dust flows. As t increases, more and more of the energy moves into the central 
region: for t > a\/2, the energy in the central region reads 
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Figure 5: Paths (|55p and of the fluid elements moving initially to the left with 2:+ = xf . In the 
Minkowski coordinates {t, x), the fluid comes momentarily to rest on the light cone of the origin, and 
it stays at rest everywhere within the interaction diamond. 



Note that as T^jy is conserved, i?ccntrai(i) can be alternatively calculated as the integral of 
coshry times the comoving energy density e(r, rf) (j53p over the constant r surface that meets 
the constant t surface on the boundary of the central region, tV2/a - 1 = TV(2a2): 



-E^central(T") = 2 / di] coshf] Te{T > aV2, ry) = 1 2" 

Jo 327rG3 ay ^ 



(58) 



The factor cosh 77 here comes from the relative velocity of the flow lines, orthogonal to the 
constant r surfaces, and the vector with respect to which the conserved energy was defined. 

Now, to what extent can we recover time-dependent thermodynamics for the fluid flow after 
the beams have collided? The energy density is given by ()52p . or equivalently (j53p . and the 
equation of state is e = p, as is expected of a (l+l)-dimensional situation. The crucial step to 
thermodynamics is to find a reasonable notion of temperature. In the fragmentation regions 
the situation is unclear, for reasons discussed in subsection [331 In the central region, however, 
the flow is Bjorken similarity flow. If boundary effects can be neglected, the arguments of 
subsection 13.11 give this flow the r-dependent temperature and entropy density 

_ ^ _ ^ ^ 

-'central — ~^ ; •^central — "TPi ' l""^/ 

ZTTT 4Gr3 r 
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Neglecting the boundary effects is likely to be justified provided the proper length tAt] of 
a constant T curve through the central region is much larger than the length scale ^^^j^trai 
associated with the temperature, and for values of \r]\ that are not close to the boundary of 
the central region. These conditions read 

^ «;log— |r/|<Clog— (60) 



or equivalently 



QgTT/VAf ^ Vx+x-, a<Cx=^. (61) 

Similarly, if the boundary effects in the interaction diamond are neglected, the arguments of 
subsection 13.31 give the static flow therein the constant temperature and entropy density 

diamond n /TT ' ■'diamond — r-ri V'-'^j 

2\'2iTa 4(_r3 2 

and neglecting the boundary effects is likely to be justified provided the diamond is much 
larger than the length scale T^J^^mond associated with the temperature and we are not close 
to the boundaries of the diamond, which conditions read 

1 < M, |x=^ - ia| < a. (63) 

The formulas (j59p and (j62p for the temperature and the entropy density are thus likely 
to be reliable whenever M ^ 1, which can be regarded as the condition that the bulk is 
semiclassical [35], and we are not close to the boundary of respectively the central region and 
the interaction diamond. Note that these conditions do not depend on how a compares with 
the other length scales of the problem, C and C,\fM. As a consistency check, we note that in 
the limit of large M we have \/M ~\ w \/M, and formulas ([59]) and ([62]) for the temperature 
and the entropy density then agree at the point (r, r\) = (a-v/2, 0), at which the central region 
and the interaction diamond meet. 

The model provides thus a consistent thermodynamical picture both in the hard core 
interaction region immediately after the collision, where the temperature and entropy density 
are approximately constant, and at late times, where the temperature falls off as l/r and 
the fluid is a perfect fluid in adiabatic expansion [9]. In Section [5] we shall compare this 
(l+l)-dimensional thermodynamics with what would be expected in 3+1 dimensions in the 
classical gluon field model with currents of nuclei as sources. 



5 Ensemble of classical gluon fields 

Contrast now the above with the modeling of heavy ion collisions using ensembles of classical 
color fields [34J. One solves the color fields A^(t,x) in suitably chosen gauges from the 
Yang-Mills equations 

[z)^,F^ = J^ (64) 

where the source current represents the two nuclei moving along the light cone, 

= 5'^+5(x-)p(i)(xt) +<^'^-5(x+)p(2)(xt). (65) 
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This current exists for all x and has the nuclear transparency built in it. The two color 
charge densities P(^jti-^{xt) are, independently for the two nuclei, drawn from a Gaussian (in 
the simplest version of the model) random ensemble, 

{pf^^{xT)pU{yT)) = 9^fiUS'^'S^{^T-yT), m = l,2, (66) 

where g'^p is a parameter describing the transverse density of color charges related [36j to the 
saturation scale by ^ 0.36(5^^)^ 0.05GeV^^^/^(\/s/Q<j)°-^. Ah information on the size 
of the nucleus and collision energy is built in Qs- The computed expectation values (A^j^A^) 
can be converted to multiplicities and initial energy densities. 

Numerical results are reported, for example, in |37j . One finds the initial energy density 
e(r, = 0) ~ for t<1/Qs', thermalisation takes place at later times. In view of the different 
dimensionalities the comparison is only qualitative, but with the identification a ~ 1/Qs the 
behavior is similar to that in Eq. (j54p . We thus suggest this interpretation for the parameter 
a introduced above. The gravity/CFT duality picture has also thermalisation built in. 

6 Conclusions 

We have in this paper presented an exact AdSs gravity solution which represents in the 
bulk a collision of two plane waves and on the conformal boundary a collision of two (1+1)- 
dimensional extended objects. We chose the incoming wave profiles to consist of a hard core 
"nucleus", of boundary length scale o that can (but need not) be comparable to the AdS 
length scale C in the bulk, and a "tail" with a 1/x^ falloff . In the region after the collision we 
then identified an interaction diamond where the hard core nuclei overlap, a central region 
where the tails overlap, and two interpolating fragmentation regions. 

We saw that the boundary matter recoils in the collision: the conformal matter is thus ex- 
tremely opaque. The interaction diamond corresponds, qualitatively, to the part of collision 
dynamics described by classical Yang-Mills field models with the indentification a ~ 1/Qs, 
where Qs is the saturation scale. There is also some qualitative similarity to Landau's hy- 
drodynamic model [291 EH] • 

We identified the bulk duals of the central region and the interaction diamond as parts of a 
spinless BTZ black hole spacetime, with the angular dimension unwrapped, and in each case 
the local invariance of the boundary flow under a spacelike Killing vector was used to identify 
a horizon and its Hawking temperature in the bulk. When the colliding matter is sufficiently 
energetic, these bulk black holes are semiclassical (M 1), and we used the bulk-boundary 
correspondence to induce from the bulk temperature a temperature to the boundary flow. 
We argued that in the limit Af S> 1 the transition effects between the regions are negligible 
provided one is not close to the boundary of the region. In this approximation the flow in 
the interaction diamond the is static with constant temperature and entropy density, and the 
flow in the central region is Bjorken similarity flow of a perfect fluid in adiabatic expansion, 
with temperature inversely proportional to the proper time 

There are clearly many open issues to be studied. Within our colliding plane wave solution 
in the AdSa bulk, one question is how to introduce thermodynamics even when the boundary 
flow is not invariant under a spacelike Killing vector. The answer would open an avenue 
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for investigating various phenomenologically-motivated thermalisation scenarios by suitable 
choices of the incoming wave profiles. The main issue of interest would clearly be to extend 
the work to more spatial dimensions, where gravitational bulk dynamics can be expected to 
give more structure to the boundary theory. 
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